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Abstract: We review recent results on new physical models constructed as VT- 

1 

symmetrical deformations or extensions of different types of integrable models. We present 
non-Hermitian versions of quantum spin chains, multi-particle systems of Calogero-Moser- 
Sutherland type and non-linear integrable field equations of Korteweg-de-Vries type. The 
Oh quantum spin chain discussed is related to the first example in the series of the non-unitary 

models of minimal conformal field theories. For the Calogero-Moser-Sutherland models 
we provide three alternative deformations: A complex extension for models related to all 
types of Coxeter/Weyl groups; models describing the evolution of poles in constrained 
real valued field equations of non-linear integrable systems and genuine deformations 
Q\ ' based on antilinearly invariant deformed root systems. Deformations of complex nonlin- 

ear integrable field equations of KdV-type are studied with regard to different kinds of 



"PT-symmetrical scenarios. A reduction to simple complex quantum mechanical models 
CD ■ currently under discussion is presented. 



1. Introduction 



Until fairly recently Q non-Hermitian systems have been mostly viewed as not self-consistent 
descriptions of dissipative systems. However, in constrast to the previous misconception 
it is by now well understood that Hamiltonians admitting an antilinear symmetry may 
be used to define consistent classical, quantum mechanical and quantum field theoretical 
systems. Various techniques have been developed to achieve this. Central to this is con- 
struction of metric operators such that certain quantities in the models can be viewed as 
physical observables |2|, [3], ||, ||, ^, [?], ||, [| 1C, 11 1. In particular, it was found that such 



type of models models possess real energy spectra in large sectors in their parameter space, 
despite being non-Hermitian. The explanation for this property can be traced back to an 



observation made by Wigner more than fifty years ago [12], who notices that operators 



invariant under antilinear transformations possess either real eigenvalues or eigenvalues 
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occurring in complex conjugate pairs depending on whether their eigenfunctions also re- 
spect this symmetry or not, respectively. A very explicit example of such a symmetry is a 
simultaneous parity transformation V and time reversal T ■ This PT-symmetry is trivially 
verified for instance for Hamiltonian operators H, but less obvious for the corresponding 
wavefunctions if) due to the fact that often they are not known explicitly. When 

[H,PT\=0 and PTV> = V> (1-1) 

hold one speaks of a PT-symmetric system, but when only the first relation holds one 
speaks of spontaneously broken P7~ -symmetry and when none of the relations in ( |1.1| ) 
holds of broken P7~-symmetry. Here we will view the PT-operator in a wider sense and 
refer to it loosely as VT even when it is not strictly a reflection in time and space, but 
when it is an antilinear involution satisfying 

VT (an/) + P<j>) = a*VT^ + fi*VT4> for a, e C, VT 2 = I. (1.2) 

Very often synonymously used, even though conceptually quite different, are the no- 
tions of quasi-Hermiticity |[3|, 14, || and pseudo-Hermiticity |[^, 16, |4|. These concepts 



refer more directly to the properties of the metric operator and their subtle difference is 
often overlooked, even though this is very important as they allow for different types of 
conclusions. In the quasi-Hermitian case the metric operator is positive and Hermitian, 
but not necessarily invertible. It was shown [13, 14, ||] that in this case the existence 



of a definite metric is guaranteed and the eigenvalues of the Hamiltonian are real. The 
pseudo-Hermitian scenario, that is dealing with an invertible Hermitian, but not necessar- 
ily positive metric, is less appealing as the eigenvalues are only guaranteed to be real but 
no definite conclusions can be reached with regard to the existence of a definite metric. 
Thus in this latter case the status and consistency of the corresponding quantum theory 
remain inconclusive. 

Even though some fundamental questions remain partially unanswered, such as the 
puzzle concerning the uniqueness of the metric or the question of what constitutes a good 
set of ingredients to formulate a consistent physical theory, the understanding is general 
in a very mature state. So far it could be used to revisit some old theories, which had 
either been discarded as being non-physical or had considerable gaps in their treatment, 
and put them on more solid ground. Another interesting possibility which had opened up 
through these studies is the formulation of entirely new models based on non-Hermitian 
Hamiltonians which however possess the desired PT-symmetry In other words, one may 
use the PT-symmetry to deform or extend previously studied models and thus obtain large 
sets of entirely unexplored theories. In principle, this kind of programme can be carried 
out in any area of physics. Here we will explore how these ideas can be used in the context 
of integrable models. We will not report here on how well established methods from 
integrable systems can be applied to study non-Hermitian quantum mechanical models 



1 17], even though we will report some scenarios in which they naturally emerge as reduced 
integrable systems |jT8| . Instead we present here how these ideas have been used so far to 
formulate and study new models previously overlooked as they would have been regarded 
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as non-physical due to their non-Hermitian nature. We present results on standard types 
of integrable models, a quantum spin-chain, multi-particle systems of Calogero type and 
nonlinear wave equations of KdV-type. 

The construction principle is fairly simple. Identifying some anti-symmetric operators 
O in the system, we seek a deformation map 5 e of the form 



VT : O h+ -O => 5 e : O ^ -i(iO) £ , (1.3) 

with e being a deformation parameter such that the non-deformed model is recovered in 
the limit e — > 1. Alternatively one can also just add 'PT-symmetric terms to the original 
system and regard them as perturbations. 



2. PT-symmetrically deformed quantum spin chains 

Quantum spin chains constitute a good starting point since, being just finite matrix models, 
they can be viewed in many ways as the easiest integrable models. We present here a model 
which has been considered first by von Gehlen @, that is an Ising quantum spin chain 
in the presence of a magnetic field in the z-direction as well as a longitudinal imaginary 
field in the x-direction. The corresponding Hamiltonian for a chain of length N acting on 
a Hilbert space of the form (C 2 )® N is given by 



1 - 

H{X,k) = ~2^K + Xa j a j+i + iK(jX j)' A,kGM. (2.1) 



We used the standard notation for the 2^ x 2 -matrices a^' y,z = I ® I ® . . . ®a x,y ' z ® . . . ( 
I ® I with Pauli matrices 






**=H' ° >= l'»J' ^l"-^' <22) 

describing spin 1/2 particles as i-th factor acting on the site i of the chain. This model 



is of interest as it can be viewed [ 20 1 as a perturbation of the A^s^-model in the Ai pq - 
series of minimal conformal field theories . It is the simplest non-unitary model in this 
infinite class of models, which are all characterized by the condition p — q > 1 and whose 
corresponding Hamiltonians are all expected to be non-Hermitian. The "P7~-symmetry of 



the model was exploited in [22]. 



2.1 Different versions of "PT-symmetry 

Let us first identify the PX-symmetry for the Hamiltonian (|2.1|). Non-Hermitian spin 



x,y,z 



chains have first been studied in this regard in 23], where the parity operator V' : o\ 

was interpreted quite literally as a reflection about the center of the chain. Viewing 
T as a standard complex conjugation V'T is then easily identified as a symmetry of the 
XXZ-spin chain Hamiltonian Hxxz- However, it is seen immediately that this operator 
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is not a symmetry of the Hamiltonian H{\, k) in (2.1). Defining instead 22] the operator 

N 



V :-- 



n 

i=i 



o-J 



with P 2 = I 



as an analogue to the parity operator, we may carry out a site-by-site reflection 



7> : (af,of,o?)^(-af,-af,of) and T : (of , o? , a* 



(2.3) 



(2.4) 



It is then easy to verify that this operator is a symmetry of H (A, k), i.e. we have [PT, H(\, k) 
0. Clearly this PT-operator acts antilinearly satisfying (|1.2| ) and is therefore a viable can- 
didate for our purposes. In analogy to ( [2 .3D it is then also suggestive to define 



N 



N 



V x :=l[af and 7>„ := IH' 



(2.5) 



i=i 



which act as 

V x/y : (of, of, of) (±of, T of,-of). (2.6) 

One can verify that [V x / y T, H(X, «)] / and [P^/yT, ffxiz] = 0. Similar properties can 
be observed for the non-Hermitian quantum spin chain [E4|] 



N 



HDG = ^2 ( K ^ a i a i+l + K x a i + K V a i) > 



(2.7) 



i=l 



with k zz G R and k v G C. Clearly when k x or ^ R the Hamiltonian Fog is 
not Hermitian, but once again one can find suitable symmetry operators. We notice that 
[P'T, H] ^ 0, whereas [PT, H] = for k x ,k v G zR and [P x / y T, H] = for kwj, G R, w^a,. G 



Below we will encounter further ambiguities in the definition of the antilinear symme- 
try, which will all manifest in the non-uniqueness of the metric operator and therefore in 
the definition of the physics described by these models. For the Hamiltonians Hxxz and 
Hog the consequences of this fact are yet to be explored. 

2.2 The two site model 

It is instructive to commence with the simplest example for which all quantities of interest 
can be computed explicitly in a very transparent way. We specify at first the length of the 
chain to be N = 2 and without loss of generality fix the boundary conditions to be periodic 

rr x 



af. The Hamiltonian fl2,l|) then acquires the simple form of a non-Hermitian 



(4 x 4)-matrix 

H 2 {\,k) 



I + I ® a z + 2Xa x <g> a x + in (I <g> a x + a x <g> I)] 



/-Iff A\ 

| A | 

%k \ n ift 

T A u T 

\ in IK 1 



(2.8) 



(2- 



- 4 - 



PT- symmetric deformations of integrable models 



The characteristic polynomial for ( |2.9[) factorizes into a first and a third order polynomial 
such that the eigenvalues acquire a simple analytic form. Defining the domain 

Up? = {A, k : k 6 + 8A 2 k 4 - 3k 4 + 16A 4 k 2 + 20AV + 3k 2 - A 2 - 1 < 0} (2.10) 

in the parameter space, the PX-symmetry is unbroken in the sense described by when 
(A, k) 6 U-p-y. The four real eigenvalues are evaluated in this case to 



e\ = A, ei = 2q2 cos (|) - |, 



e 3/4 = 2g2 cos (f + vr =F f) 



A 
3 > 



(2.11) 



where 



9 = arccos 



7 3/2 



- (3 + 4A 2 - 3k 2 ) , r = A(l8^ + 8A 2 + 9). (2.12) 
9 27 



We depict the eigenvalues in figure |T] for some fixed A or k and varying k or A, respectively. 




0.0 0.1 



0.0 0.1 



0.3 0.4 



Figure 1: Unavoided level crossing: eigenvalues as functions of A (k) for fixed k (A) for i?2(A, n). 



We observe the typical behaviour for "PT-symmetric systems, namely that two eigen- 
values start to coincide at the exceptional point |25j when k and A are situated on the 
boundary of U-p-y. Going beyond those values, the "PT-symmetry is spontaneously broken 
and the two merged eigenvalues develop into a complex conjugate pair. This is of course 
a phenomenon prohibited for standard Hermitian systems by the Wigner-von Neumann 
non-crossing rule |]26| . 

For the Hamiltonian ( |2.8| ) one can compute explicitly the left |<3? n ) and right eigenvec- 
tors l^n); forming a biorthonormal basis 



and V |*„)(*„| 

*n 



I 



(2.13) 



and verify that indeed for the spontaneously broken regime the second relation in ( |1.1| ) 
does not hold, see |^2[ for the concrete expresssions. We have then all the ingredients to 
compute the metric operator p and define the inner product (.\.) p ■= {-\p-) with regard to 
which the Hamiltonian (|2.8|) is Hermitian 



(2.14) 
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Computing the signature s = (si, s%, . . . , s n ) from 

V\<5> n ) = S n \^ n ) Withs n = ±l (2.15) 

we may evaluate the so-called C-operator introduced in 

C:=V a n |* n )(* n |, (2.16) 
* — 'n 

and hence the metric operator p, which also relates the Hamiltonian to its conjugate 

p:=VC, H^p = pH. (2.17) 



The explicit expressions can be found in [22], from which one can verify explicitly that 
the metric operator is Hermitian, positive and invertible. Thus the Hamiltonian (2.8) is 



quasi-Hermitian as well as pseudo-Hermitian. From the expression for p we can also obtain 
the so-called Dyson map rj [27], by taking the positive square root r\ = yfp. This operator 
serves to construct an isospectral Hermitian counterpart h to H by its adjoint action. For 
fO|) we find 

h 2 (\, k) = r]H 2 (\, k)!]' 1 = ^ v s o- s ®o- s + p z (o- z ®l + l®o- z ). (2.18) 

s=x,y,z 

The constants u x , u y , v z and p z can be found in [^]. 

2.3 Perturbative computation for the N-site model 

It is clear that when proceeding to longer spin chains it becomes increasingly complex to 
compute the above mentioned operators such that exact computation become less transpar- 
ent and can only be carried out with great effort. However, we may also gain considerable 
insight by resorting to a perturbative analysis. For this purpose we separate the Hamilto- 
nian into its Hermitian and non-Hermitian part as H(X,k) = ho(X) + inhi, where ho and 
hi are both Hermitian, with k being a real coupling constant as for instance introduced in 
(|2.1j) . Assuming next that the inverse of the metric exists and that it can be parameterized 
as p = e q , the second equation in ( p,17|) can be written as 

& = e<>He-i =H+[q,H) + ~[q, [q, H}] + ^[q, [q, [q, H]}] + ■■■ (2.19) 
Presuming further that the metric can be perturbatively expanded as 

oo 

g = Xy fc ~Wl. (2.20) 

k=l 

we obtain the following equations order by order in k 

[h , qi ] = 2ih x , (2.21) 

= ^[<7l,[<7iA]], (2-22) 

It % 

[V*] = g[gi, [Q3,hi]] + ~[q 3 , [qiMW ~ ggQ^i. [ffii bi,fri]]]]- (2-23) 
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It is clear from ( 2.21| )-( |2.23 ) that at each order one new unknown quantity enters the 



computation for which we can solve our equations, i.e. in (2.21) we solve for q\ for known 
ho and hi, in ( |2.22| ) for (73, in ( |2.23| ) for q$, etc. This process can be continued up to any 
desired order of precision, see || for further general details on perturbation theory. 

Proceeding in this manner we compute the Dyson operator as described above and 
determine the Hermitian counterpart. For iV = 3we obtain 

u — ,.3 o3 _|_, 3 c-3 , 3 c-3 , 3 o3 , 3 o3 , 3 o3 , 3 c-3 (9 9A\ 

11 h L xx u xx ~ H'yy u yy ' H'zz u zz ~ h L z u z ~ H'xxz u xxz ' h L yyz u yyz ~ H'zzz kJ zzz> V^-^l 

where for convenience we introduced a new notation 

N 

S a 1 a 2 ...a p -=^°t 1 ° a k\i---e a k +P -v for a,i = x,y, z,u; i = 1, . . . ,p < N. (2.25) 

fc=l 

The coefficients p xx , . . . , \i zzz are real functions of the couplings A and k. We denote here 
a u = I to allow for the possibility of non-local, i.e. not nearest neighbour, interactions. In 
fact, they do occur when we increase the length of the chain by one site. For N = 4 we 
compute 

U _ .,4 c4 , 4 c4 i 4 rr4 , 4 c4 , 4 r-4 , 4 c4 , 4^4 

rn'-'ii ' xx^xux ' t^yy yy yy yuy ' H'zzuzz ' zz^zuz ' t^z^z 

_)_ i y4 (oA , c4 \ , ..4 c-4 , 4 ( a A , c4 \ . ..4 c-4 

T rrazV u M2 ^ ^zxxl ' r^xzx^xzx ' r^yyzx^yyz ' ^zyyJ ' H'yzy^yzy 

+ „ 4 Q 4 + ,,4 C-4 , 4 c4 , 4 c-4 , 4 r-4 

' ^222^22:2 ^ H'xxxx'-'xxxx ' r^yyyy yyyy ' H'zzzz'~ J zzzz < r'xxyy'-'xxyy 

_l„ 4 Q 4 _l ,, 4 Q4 _j_ ,,4 c-4 , 4 c-4 , 4 o4 (9 96^ 

< r^xyxy^xyxy ' zzyy^ zzyy ' r b zyzy u zyzy ' r^xxzz^xxzz ' r*'xzxz'- J xzxz' 

We observe that the first non-local interaction terms proportional to S xux , Sy Uy and S zuz 
emerge in this model. Thus we encounter a very typical feature of non-Hermitian VT- 
symmetric Hamiltonian systems, whereas the non-Hermitian Hamiltonian is fairly simple 
its Hermitian isospectral counterpart is quite complicated involving non-nearest neighbour 
interactions. An additional feature not present for chains of smaller length is the fact that 
some of the A-dependence of the coefficients f/, xx , . . . , pt xzxz is no longer polynomial and 
gives rise to singularities. 

This models exhibits the basic feature, but clearly there is plenty of scope left for 
further analysis. More explicit analytic formulae should be computed for rj, p and h for 
longer chains, models with higher spin values should be considered and further members of 
the class belonging to the perturbed .M Pi <j-series of minimal conformal field theories should 
be studied. Interesting recent results on other non-Hermitian quantum spin chains may be 
found in [M|l. 



3. PT-symmetrically deformed Calogero type models 

'PT-deformed versions of multi-particle systems of Calogero type have been obtained so 
far in three quite different ways, as simple extensions, as constrained field equations or as 
genuine deformations. 
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3.1 Extended Calogero-Moser-Sutherland models 

The most direct and simplest way to obtain PT-symmetrically extended versions of Calogero- 
Moser Sutherland models is to add a PT-symmetric term to the original model as proposed 



in |0] 



UvTCMS = \p 2 + 9a v ( a ■ l) + ^9af(a ■ q)a ■ p, (3.1) 

with coupling constants canonical variables q,p G for an (£ + l)-dimensional 

representation of the roots a of some arbitrary root system A, which is left invariant under 
the Coxeter group. The potential may take on different forms V(x) = f 2 (x) defined by 
means of the function f(x) = 1/x, f(x) = 1/sinhx or f(x) = 1/ sinx. The model in ( p.l| ) 
is a generalization of an extension of the A(_ and -B^-Calogero model, i.e. f(x) = 1/x, for 
a specific representation of the roots as suggested in ||3l"l , [32] , |33|| . The PT-symmetry of 



rl-PTCMS is easily verified. In [ 30 1 it was shown that for f{x) = 1/x one may re-write the 
Hamiltonian in (|3.1[) , such that it becomes the standard Hermitian Calogero Hamiltonian 
with shifted momenta p — > p + i/i 

HvTCMS = \(P + + \Y. 9l V ( a ■ (3-2) 

aGA 

where [i := 1/2 J2 a eA 9af( a ' Q) a an d the coupling constants have been redefined to := 
gl + ct 2 s gl for a £ A s and g^ := gf + afgf for aGA/, where A; and A s refer to the root 
system of the long and short roots, respectively. This manipulation is based on the not 
obvious identity fj? = a 2 s g 2 s X^aeA s V( a ' l) + a fdf S a eA ; V( a " q)> which is only valid for 
rational potentials. Even then it has not been proven yet in a case independent manner, 
but verified for many examples on a case- by-case basis [3C]. 



For the rational potential it is straightforward to obtain the Dyson map rj = e 



which relates the standard Hermitian Calogero model to the non-Hermitian model (3.1) by 
an adjoint action T-L-ptc = "ri~ l 7~LcV- The integrability of the rational version of 'H.vtcms 
follows then from the existence of the Lax pair Lc and Mq obeying the Lax equation 
L-pTC = [L-pTC 'j M-pfc], which maybe obtained from the standard Calogero Lax pair ]34| 
as L VTC {p) = rf x L c (j))ri = L c {p + iy) and M VTC {p) = ■q~ 1 M c {p) , q = M c {p + in). 
Expanding the shifted kinetic term in (|3.2| ) we obtain 

KvTCMS = \p 2 + 9l v ( a ■q) + ifJ>-p - ^/i 2 - (3.3) 

aeA 

By the reasoning provided, it follows that this model is integrable for all of the above stated 
potential, whereas the model without the /i 2 -term is only integrable for rational potentials. 

3.2 From constraint field equations to "PT-deformed Calogero models 

Another more surprising way to obtain particle systems of complex Calogero type arises 
from considering real valued field solutions for some nonlinear equations. Making an Ansatz 
in form of a real valued field 

^"-iS^"^)' A6K ' <3 ' 4) 
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it was shown more than thirty years ago 36] that this constitutes an ^-soliton solution 
for the Benjamin-Ono equation 

u t + uu x + \Hu xx = 0, (3.5) 
with Hu(x) denoting the Hilbert transform Hu{x) = — j^dz, provided the poles z k 



in (3^) obey the complex ^-Calogero equation of motion 

.. A 2 



Zk = Y ^^ Zj ~ Zk) 3 ' Zk e C ' 



Clearly for different types of nonlinear equations the constraining equation might be of a 
more complicated form. However, we may consistently impose additional constraints by 
making use of the following theorem of Airault, McKean and Moser [37|: 
Given a Hamiltonian H(x%, . . . , x n , x%, . . . , x n ) with flow 

Xi = dH/dii and Xi = —dH/dxi i = l, ...,n (3-7) 

and conserved charges Ij in involution with H, i.e. vanishing Poisson brackets {Ij,H} = 0. 
Then the locus of grad I = is invariant with regard to time evolution. Thus it is permitted 
to restrict the flow to that locus provided it is not empty. 
Making now the Ansatz 

t 

v{x,t) = X^2 (x - z k (t)y 2 , AeK (3.8) 

k=l 

one can show that this solves the Boussinesq equation 

v tt = a(v 2 ) xx + bv xxxx + v xx a,kl. (3.9) 
if and only if b = 1/12, A = —a/2 and the poles z k obeys the constraining equations 

z k = 2j2(zj ~ ^k)' 3 Zk = -^f^; ( 3 - 10 ) 

z 2 k = l- ~ Z *Y 2 S rad ( J 3 - h) = 0. (3.11) 

Here I3 = Ylj=i[Zj + Tl/k^j Zj{ z j ~ z k) 2 ] and I\ = J2 £ j=i z j are two conserved charges 
in the ^-Calogero model. Thus in comparison with the previous example ( |3 . 4| ) - ( pT5| ) we 
have to satisfy an additional constraints ( |3.11| ) besides the equations of motion of the Ag- 
Calogero model. However, according to the above theorem this is still a consistent system 
of equations, provided the equations ( 3.10 ) and ([3.11 ) possess any non-trivial solution. 



Only very few solutions have been found so far. The simplest two-particle solution was 
already reported in [j37|l 



zi = k + y/(t + k) 2 + 1/4, z 2 = k - ^{t + Rf + 1/A. (3.12) 
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In this case the Boussinesq solution acquires the form 

v(x t) - 2A {X ~ K)2 + {t+ ~ K? + 1/4 (3 13) 

Note that v(x, t) is still a real solution. However, without any complication we may change 
k and k to be purely imaginary in which case, and only in this case, ( |3. 13 ) becomes a 
solution for the "PT-symmetric equation (3J3) in the sense that VT : x — > —x,t — > —t and 
v — > v. A three particle solution was reported in f}q |, which exhibits an interesting solitonic 
behaviour in the complex plane. In that case no real solution could be found and once 
again one was forced to consider complex particle systems. For more particles, different 
types of algebras or other types of nonlinear equations these investigations have not been 
carried out yet. 

3.3 Deformed Calogero-Moser-Sutherland models 

Let us now consider the CMS- models with an additional confining potential 

K'vTCMS = Y + l6 ^ ( a ' ( l) 2 + lY,9 a V(.a-q), m,g a eR, (3.14) 

aeA s aeA 

and also deforme the coordinates q —> q. Considering at first the ^2-case for a standard 
three dimensional representation for the simple ^-roots ct\ = {1, —1,0}, a 2 = {0, 1, — 1}, 
we deformed the coordinates as 

Qi — > Qi = Qi coshe + iV3(q 2 — q 3 ) sinhe (3.15) 
<?2 — > q~2 = <?2 cosh e + iV3(q 3 - qi) sinhe (3.16) 
<?3 — > h = ^3 coshe + iV3(q\ - 92) sinhe (3-17) 

such that the relevant terms in the potential become 

OL\-q = qn coshe i=(qi3 + Q23) sinhe, (3.18) 

V3 

«2 • q = <?23 coshe 7= (521 + <?3i) sinhe, (3.19) 

v3 

(ai + a 2 ) • q = Qi3 coshe + —^(q%2 + 932) sinhe, (3.20) 

v3 

with the abbreviation qij := qi — qj. We observe for this example the following antilinear 
involutory symmetries 

51 : gi« g 2 , <?3 <?3, i -> (3.21) 

5 2 ■■ q2 -H- q 3 , qi -H- q x , i ->• -%. (3.22) 



At this stage this deformation appears to be somewhat ad hoc. In fact, it arose [39, 40] 
from the physical motivation to eliminate singularities in the potential when solving the 
separable Schrodinger equation for the Hamiltonian W-pfcMS- ^ was n °ted that the new 
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non-Hermitian model could be defined on less separated configuration space. Whereas in 
general one had to restrict the models to distinct Weyl chambers and analytically continue 
the wavefunctions across their boundaries with the inclusion of some chosen phase, this 
is no longer necessary in the deformed models. In addition, the new models possess a 
modified energy spectrum with real eigenvalues, which we attribute to the fact that the 
theory is invariant with respect to the antilinear transformations S\ and 52- Motivated by 
this success one may attempt to find a more direct systematic mathematical procedure to 
deform the coordinates rather than the indirect implication resulting from the separability 
of the Schrodinger equation. In any case, the latter approach would be entirely unpractical 
for models related to higher rank Lie algebras. 

We notice first that the Hamiltonian ( |3. 14 ) also results from deforming the roots 



involved. For the A2-case we may take the simple roots 



1 

dii = ai coshe + i— — sinhe(ai + 2a2), (3.23) 
v3 

a>2 = ct2 coshe — i—j= sinhe(2o;i + 02). (3.24) 
v3 



and re- write ( 3.14 ) equivalently as 

2 ^2 ^ 
n'-prcMS = y + XJ (« • 1? + 5 S 9&V ^ ' m ' 9& G M " (3 ' 25) 

Now the symmetries fl3,21| )-( |3,22| ) can be identified equivalently for the roots. We note 

<rf : &x <-> — 5i, a 2 O ai + a 2 O qi O q 2 , q3 «-> q3, i ~ > ~h (3.26) 
o\ : «2 -«2> oil "i + "2 ^ 92 93, 91 ^ qi, i -> —i. (3.27) 

This observation has been taken as the basis for the formulation of a systematic construction 
procedure leading to antilinerly invariant, and therefore potentially physical, models [41 



42, ffij. The dynamical variables, or possibly more general fields, appear in the dual space 
of some roots with respect to the standard inner product. Since these root spaces are 
naturally equipped with various symmetries due to the fact that by construction they 
remain invariant under the action of the entire Weyl group W, it is by far easier and 
systematic to identify the antilinear symmetries directly in the root spaces rather than in 
the configuration space. Once they have been identified they can be transformed to the 
latter. 

The aim is therefore to construct complex extended antilinearly invariant root systems 
which we denote by A(e). The proposed procedure consists of constructing two maps, 
which may be obtained in any order. In one step we extend the representation space A of 
the standard roots a from W 1 to C n = M n © iW 1 . This means we are seeking a map 

6 : A -> A(e), a i-> a = 6 £ a, (3.28) 

where a = {on, . . . , o>g}, A C W l , A(e) C C n and n is greater or equal to the rank i 
of the Weyl group W. The complex deformation matrix 9 £ introduced in ( [3.28D depends 
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on the deformation parameter e in such a way that lim^o 9 £ = I. The deformation 
is constructed to facilitate the root space A with the crucial property for our purposes, 
namely to guarantee that it is left invariant under an antilinear involutory map 

w : A(e) A(e), a ^ u&. (3.29) 

This means the map in ( p.29| ) satisfies w : a = + /U2«2 ^ M* WQ( i + A t 2 u;Q! 2 for /ii, 

/i2 £ C and also o = I. In order to facilitate the construction we make the further 
additional assumptions: 

(i) The operator u decomposes as 

ui = tuj = Cut, (3.30) 

with oj G W, cD 2 = I and r being a complex conjugation. This will guarantee that w 
is antilinear. 

(ii) There are at least two different maps cjj with i = 1, ...,« > 2. This assumption 
simplifies the solution procedure. 

(iii) There exists a similarity transformation of the form 

u>i := 9 £ Cui9~ 1 = T(2)i, for i = 1, . . . , k > 2. (3.31) 

(iv) The operator # £ is an isometry for the inner products on A(e), such that 

1 and det0 e = ±l. (3.32) 



This assumption is motivated by the desire to keep the kinetic term of the Calogero 
model undeformed. 

(v) In the limit e->0we recover the undeformed case 

lim 6 £ = L (3.33) 

Clearly one could modify or entirely relax some of the constraints (i)-(v), e.g. it might 
not be desirable in some physical application to preserve the inner products etc. However, 
it turns out that this set of constraints is restrictive enough to allow for the construction 
of solutions for £ with only very few free parameters left. 

With our applications to physical models in mind, i.e. exploiting here the equivalence 
of ( |3.14| ) and ( p. 25 ), we would also like to construct a dual map 5* for 5 acting on the 



coordinate space with q = {qi, . . . , q n } or possibly fields. We therefore define 

5* : R n ^ A*(e) =M n eilfP, x^x = 6*x, (3.34) 

denoting quantities in and acting on the dual space by *. Thus assuming 9 £ has been 
constructed from the constraints (i)-(v), we may obtain 6* by solving the £ equations 

{an ■ x) = {{9 £ a) i ■ x) = {oti ■ 6 £ x) = (aj • x), for i = 1, . . . ,£, (3.35) 
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involving the standard inner product. This means (9*) l oti = (9 £ a) i . Note that in general 
0* ^9*. Naturally we can also identify an antilinear involutory map 

w* : A*(e) -»■ A*(e), x ■-> (3.36) 

corresponding to w but acting in the dual space. Concretely we need to solve for this the 
K x £ relations 

(u)i5t) ■ ■ x = ctj • co*x, for i = 1, . . . k; j = 1, . . . ,£, (3.37) 
for u* with given o;j. 



In [41, 42, many solutions to the set of constraints (i)-(v) were constructed. A 
particular systematic construction can be found when we take k = 2 in the requirement 
(ii) and identify oj\ = c_ and ui2 = cr + . The maps a± factorize the Coxeter element in a 
unique way 

a := a-cr + , with a± := j^J cij, (3.38) 

i€V± 

where the Oi are simple Weyl reflections Oi{x) := x — 2(x-oti/af)ai associated to each simple 
root for 1 < i < t = rankW. The two sets V± are defined by means of a bi-colouration 
of the Dynkin diagrams consisting of associating values q = ±1 to its vertices in such a 
way that no two vertices with the same values are linked together. The consequence of 
this labeling is that [crj,c7j] = for i,j G V+ or i,j £ V~ such that the factorization in 



( 3.38 ) becomes unique. Clearly o± = I as required by for our construction. An immediate 



± 

consequence of (iii) is that a and 8 e commute, such that the following Ansatz captures all 
possible cases based on the assumption stated before ( |3.38 ) 



£ = y^c k (e)a k , with lim c fc (e) = <^ Cfc ( e ) G C. (3.39) 

fc=o L r 

The upper limit in the sum results from the fact that a h = I, with h denoting the Coxeter 
number. Invoking also the remaining constraints allows to determine the functions Cfc(e). 
For the ^3-Weyl group invariant root system this yields for instance the following three 
deformed simple roots 



ajeosheai + (coshe — 1 ) 0:3 — i V2V cosh e sinh y-j (ai+2a2+«3) , (3.40) 



a^2 cosh e — 1)(X2 + 2iV2V cosh e sinh y-j (a\ + 02 + 03) ; (3-41) 



ageoshea3 + (coshe — l)a\ — iv^Vcoshesinh y-j (01+202+03). (3.42) 

In some cases we were even able to provide closed formulae for entire subseries. For instance 
for A^n-i we found a closed expression for the deformation matrix in the form 



6 £ = r I + (1 - r )a 2n + i^r 2 - r {a n - a' n ) . (3.43) 

A possible choice for the function ro is ro = coshe. It was also shown in [41] that it is 
impossible to construct solutions for (i)-(v) for certain Weyl groups based on the factor- 
ization ( |3.38 ), such as for instance i?2n+i- However, in (4^, Q it was demonstrated that 
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one can slightly alter the procedure by chosing different factors instead and constructing 
solutions based on an Ansatz similar to (3.39). For B2 n +i we found closed expressions of 
the form 



OL2j-i = cosheQ2j-i + % sinhe a^'-l + 2 I for j = 1, . . . , n, (3.44) 



k=2j 



2j+2 



ct2j = cosheQ2j — i sinhe «fc + 2 2ak I for j = 1, . . . , n — 1, 

\fc=2j fe=2i+3 / 

ot£_i = coshe(o^_i + ag) — a.£ — i sinhe (a^_2 + + a^) , 



(3.45) 



In this case the dual deformation matrix which acts on the coordinates Q3.35 ) takes on a 
very familiar form and turns out to be composed of pairwise complex rotations 



(R \ 
R 



R 



\ 



with R 



cosh e i sinh e 
-i sinhe coshe 



(3.46) 



Having constructed various deformed root spaces which by construction are equipped 
with an antilinear involutory symmetry, we may then consider various models formulated in 



terms of roots, such as H'-pn-CMS defined in ( 3.25 ). We then encounter several interesting 
new features in these models. Since many of the key identities needed for the solution 
procedure are identical in terms of roots or deformed roots, we may adopt similar solution 
techniques as in the undeformed case, such separating variables. As a crucial new feature 
we find that the energy spectrum is modified and admits new real solutions when compared 
to the undeformed model. For instance, for the GVcase the undeformed energies E n = 
2 \u\ (2n + A + 1) with n G N and A G R + become 



E^m = 2\oj\ [2n + Q(nf + + m) + l] for n,m £ 



(3.47) 



with K f fl = (1 ± x /l + 4 5fl/t )/4. 

For the i?2n+i-case we can support these observations with the explicit construction of 
Dyson maps as introduced in ( |2.1§| ) and the metric operators ( 2.17| ). For the models based 
on the deformed roots ( 3.44| )-( |3T45|) the Dyson map is simply r\ = 7712^34 . . . T)u-2)U~\) with 
ijij = e - £ ( x iPj- x jPi) ; such that the metric operator becomes p = r/ 2 . 

A further novelty in the deformed models is that the wavefunctions are regularized by 
means of the deformation such that many singularities disappear. In particular this means 
that these models can be defined usually on the entire space C n , whereas the undeformed 
models could only be defined in certain Weyl chambers. The continuation from a chamber 
to its neighbouring one was achieved by introducing a phase factor by hand, thus selecting 
a particular statistics. The deformed models on the other hand have these phase factors 
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already built in as a property of the model. For instance in the ^-case we find that the 
four-particle wavefunction obeys 

^(?i,92,93,?4) = e ms ip(q 2 ,q4,(}i,q3)- (3.48) 

These properties can be read off easily from the action of the generalized 'P7~-symmetry 
on the deformed roots, translated to the dual space that is to the coordinates and then to 
the parameterization of the wavefunction. We note that the phase factor emerges as an 
intrinsic property rather than as an imposition. We illustrate the relation ( |3.48D as follows 

w x y z w x y z 

• • • • = e l7TS • • • • 

Ql 92 93 94 92 94 qi 93 

In the deformed models we can even allow some of the particles to occupy the same place 
and scatter them with single particles 

x y z x y z 

• m • = e lns • m • 

qi 92 = 93 94 92 91 = 94 93 

We can also scatter pairs of particles resulting in the exchange of one of the particles 

x y x y 

91 = 92 93 = 94 91 = 93 92 = 94 

and we may even scatter triplets with a single particle 

x y x y 

9i = 92 = 93 94 94 9i = 92 = 93 

It is clear that these models have very interesting new properties and there are still 
various open issues left worthwhile to investigate. More explicit solutions for spectra and 
wavefunctions should be constructed; the important questions of whether the deformed 
models are still integrable should be settled; Dyson maps, the metric operators and Her- 
mitian counterparts should be constructed such that more observables of the models can 
be studied. The constructed root systems could be used to formulate entirely new models 
of different kind than Calogero systems. 

4. PT-symmetrically deformed nonlinear wave equations 

The prototype integrable system of nonlinear wave type is the Korteweg-deVries equation 

m 

u t + Puu x + ju xxx = 0, /?,7 € C, (4.1) 
resulting from a Hamiltonian density 

WkivM = + (4.2) 

The system admits two different types of PT-symmetries 

VT± : x i-)- -x, 1 1-)- —t, i i-)- —i, u (->• ±u for/3,7GlR, (4.3) 
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which have only been exploited recently in 46, 47, 18 1. According to the deformation 
prescription (1.3) we can now deform the Hamiltonian density (|4.2| ) in two alternative ways 



5 £ : u x i-)- u a 



-i(iu x ) e 



or 5 £ : u ^ u £ 



-i(iu) e 



(4.4) 



respectively, depending on whether we assume u(x, t) to be PT-symmetric or PT-anti- 
symmetric. The deformed models, suitably normalised, are then defined by the densities 



n; 



7 



x ) £+ \ and HI 



6 1 + e 
with corresponding equations of motions 

u t + (3uu x + j(u Xj£ ) xx = 0, and 



(l + e)(2 + e) {m) + 2 % ' 



«t + if3u £ u x + 7n 

.7 



0. 



(4.5) 



(4.6) 



The PT-symmetry can be exploited to ensure the reality for expressions such as the energy 
on a certain interval \—a, a] 



r a r du 

E= n[u(x)]dx= <b H[u{x)}—. 

J -a JT u x 



(4.7) 



One would expect this expression to be real for the unbroken symmetric regime. However, in 
fl8[j also some unexpected cases with real energies were found for which the "P7~-symmetry is 



entirely broken (4.3), i.e. for the Hamiltonian and for its solutions. This possibly indicates 
the existence of a different realisation for the PT-symmetry operator. 





Figure 2: The two different Riemann sheets for rational solutions with broken 'PT-symmetry of 
the 'H^Yg-model with different values of purely complex initial conditions whose imaginary part is 
indicated in the figure The asymptotic conditions are chosen to be limf_ i .± 00 it(C) = (1 — *)/4, the 
speed of the wave as c = 1 and the constants of the model as j3 — 2 + 2i and 7 = 3 

A further characteristic feature of the deformed models is that, in general, one has 
to view them on various Riemann sheets. An example for a traveling wave solution pa- 
rameterized by C = x — ct, with c denoting the wave speed, for broken PT-symmetry is 
presented in figure |2|. 
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The branch cut at — oo — i/4 to (1 — i)/4 is passed from above in panel (a) to below 
in panel (b). The trajectories for the PT-symmetric and broken PT-symmetric case look 
qualitatively very similar, the major difference being that the fixed point has moved away 
from the real axis, thus leading to a loss of the symmetry. 

Viewing the systems as two dimensional models, the nature of the fixed points has 
been investigated systematically by exploiting the fact that their characteristic behaviour 
is completely classified in dependence on the different types of eigenvalues for the Jacobian. 
In |l8| it was found that they may even undergo Hopf bifurcations in these systems, passing 
form a star node over a centre to a focus. This feature was derived for the PT-symmetric 
as well as for the broken "PT-symmetric regime. In particular this also means that we 
encounter closed trajectories despite the fact that the PT-symmetry is broken. We depict 
an examples in figure for different values of c, 3 and 7. 






Figure 3: Broken PT-symmetric solution for H 4 : (a) Star node at the origin for c = 1, (3 = 2 + i3, 
7 = 1; (b) centre at the origin for c = 1, /3 = 2 + i3, 7 = —1; (c) focus at the origin for c = 1, /3 = 2, 
7 = 1-M3. 

An interesting relation between these type of deformations and some simple complex 



quantum mechanical models was pointed out in [18]. As a special case of this general 
observation we consider here the model %<^-model with Hamiltonian density 



8 7 
K2H = —u 4 + -u 2 x , 



(4.8) 



and make contact with the model studied in [48]. As explained in [ 18 [ , integrating (^ 
twice with respect to £ we obtain 



o2/ C 9 n I 4 

U C = - K2 + Kill + -U + p — U 

s 7 V 2 12 



(4.9) 



with integration constants Ki, K2 £ C. Identifying u — > x and (" — > t for the traveling wave 
equation together with the constraints 



«u = -7 r > K 2 = jEx, (3 = -375 and c = 
converts the derivative of equation (|4.9|) into Newton's equation 



-70; 



(4.10) 



9 q 

x + T + ux + gx =0, 



(4.11) 
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for the quartic harmonic oscillator of the form 



tfquartic = E X = ^p 2 + TX + yl 2 + |x 4 . (4.12) 



One may now directly translate some of the properties of the system ( |4.8| ) to the quantum 
mechanical model ( |4.12j ). The special choice K\ = r = for the integrations constants 
imply that one is considering asymptotically vanishing waves with lim^oo u(Q = and 
with Neumann boundary condition lim^oo u x (£) = V^E X where -ff qU artic = E x . Accord- 
ingly, the energy in the classical analogue of a complex classical particle corresponds to an 
integration constant in the nonlinear wave equation context multiplied by one of the cou- 
pling constants in the latter model. This is of course different from the energy as defined 
in (fh^]), which also leads to different conclusions regarding the reality of these quantities 
resulting from the various PT-symmetric scenarios. 

In a similar way, the complex seminal Q cubic harmonic oscillator 

1 1 

#cubic = 2"P 2 + ^x 2 + igx 3 , (4.13) 

treated also in |4£| simply results from the integrating the KdV-equation twice with the 
identification k\ = 0, K2 = ^E x , (3 = —i6cg and 7 = — c. 

It appears to be unlikely that the models are still integrable as in general they do not 
pass the Painleve test |30], |5lJ| . Similar studies have also be carried out for other types of 
nonlinear wave equations as for instance for deformed Ito systems in . It was even shown 
that one can T^T-symmetrically deform the supersymmetric version of the KdV-equation 



(4.1) while still preserving its supersymmetry [47|. 

Evidently many features remain still unexplored and it would be very interesting to 
extend these studies to a larger range of values for the deformation parameter, to other 
nonlinear field equations such as Burgers, Bussinesque, KP, generalized shallow water equa- 
tions, extended KdV equations with compacton solution, etc. 
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